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Abstract 

The Lagrangian of the time fractional KdV equation is derived in 
similar form to the Lagrangian of the regular KdV equation. The vari- 
ation of the functional of this Lagrangian leads to the Euler-Lagrange 
equation that leads to the time fractional KdV equation. The Riemann- 
LiouvuUe definition of the fractional derivative is used to describe the time 
fractional operator in the functional of the Euler-Lagrange formula. The 
time-fractional term of the derived KdV equation is represented as a Riesz 
fractional derivative. The variational-iteration method given by He [31] 
is used to solve the derived time- fractional KdV equation. The calcula- 
tions of the solution with initial condition Aq sech^ (ex) are carried out and 
demonstrated in 3-dimensions and 2-dimensions figures. We give the com- 
parison and estimates of the role of fractional derivative to the nonlinear 
and dispersion terms in the fractional KdV equation for unit amplitude 
Aq — 1. It worth mentioned that, exploitation of the fractional calcu- 
lus for fractional equation provides not only new types of mathematical 
construction, but also new physical features of the described phenomena. 
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1 Introduction 

Because most classical processes observed in the physical world are nonconser- 
vative, it is important to be able to apply the power of variational methods to 
such cases. A method [1] used a Lagrangian that leads to an Euler-Lagrange 
equation that is, in some sense, equivalent to the desired equation of motion. 
Hamilton's equations are derived from the Lagrangian and are equivalent to the 
Euler-Lagrange equation. If a Lagrangian is constructed using noninteger-order 
derivatives, then the resulting equation of motion can be nonconservative. It 
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was shown that such fractional derivatives in the Lagrangian describe noncon- 
servative forces [2, 3]. Further study of the fractional Euler-Lagrange can be 
found in the work of Agrawal [4-6], Baleanu and coworkers [7-9] and Tarasov 
and Zaslavsky [10, 11]. 

Recently, fractional calculus has been applied to almost every field of science, 
engineering and mathematics. The awareness of the importance of this type 
of equation has grown continuously in last decade include for visco-elasticity 
and rheology, image processing, mechanics, mechatronics, physics, and control 
theory, see for instance [12]. 

On the other hand, the Korteweg-de Vries (KdV) equation has been foimd 
to be involved in a wide range of physics phenomena as a model for the evo- 
lution and interaction of nonlinear waves. It was first derived as an evolution 
equation that governing a one dimensional, small amplitude, long surface grav- 
ity waves propagating in a shallow channel of water [13]. Subsequently the 
KdV equation has arisen in a number of other physical contexts as collision- 
free hydro-magnetic waves, stratified internal waves, ion-acoustic waves, plasma 
physics, lattice dynamics, etc [14]. Certain theoretical physics phenomena in 
the quantum mechanics domain are explained by means of a KdV model. It is 
used in fluid dynamics, aerodynamics, and continuum mechanics as a model for 
shock wave formation, solitons, turbulence, boundary layer behavior, and mass 
transport. All of the physical phenomena may be considered as nonconservative, 
so they can be described using fractional differential equations. Therefore, in 
this paper, our motive is to formulate a time-fractional KdV equation version 
using the Euler-Lagrange equation via what is called variational method [4-6, 
15]. 

Several methods have been used to solve fractional differential equations 
such as: the Laplace transformation method [16, 17], the Fourier transforma- 
tion method [16, 17], the iteration method [18] and the operational method 
[19]. However, most of these methods are suitable for special types of fractional 
differential equations, mainly the linear with constant coefficients. Recently, 
there are some papers deal with the existence and multiplicity of sohition of 
nonlinear fractional differential equation by the use of techniques of nonlinear 
analysis (fixed-point theorems, Leray-Shauder theory, Adomian decomposition 
method, variational-iteration method, etc.), see [20-24]. In this paper, the resul- 
tant fractional KdV equation will be solved using a variational-iteration method 
(VIM) [25-27]. In addition, we give the comparison and estimates of the role 
of fractional derivative to the nonlinear and dispersion terms in the fractional 
KdV equation for unit amplitude . 

This paper is organized as follows: Section 2 is devoted to describe the 
formulation of the time-fractional KdV (FKdV) equation using the variational 
Euler-Lagrange method. In section 3, the resultant time-FKdV equation is 
solved approximately using VIM. Section 4 contains the results of calculations 
and discussion of these results. 
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2 The time-fractional KdV equation 



The regular KdV equation in (1+1) dimensions is given by [13] 

^u{x,t) + Au{x,t)^u{x,t) + B ■^u{x,t)=0, (1) 

where u{x,t) is a field variable, x G R is a space coordinate in the propagation 
direction of the field and t GT {= [0,To]) is the time variable and A and B are 
known coefficients. 

Using a potential function w(a;, t), where M(a;,t) = Vx{x,t), gives the potential 
equation of the regular KdV equation (1) in the form 

Vxt {X, t) + AVx {X, t)Vxx {x,t)+B {x,t)=0, (2) 

where the subscripts denote the partial differentiation of the function with re- 
spect to the parameter. The Lagrangian of this regular KdV equation (1) can 
be defined using the semi-inverse method [28, 29] as follows. 

The functional of the potential equation (2) can be represented by 



Jiv) = j dx j dt{v{x,t)[ciVxt{x,t) + C2AVx{x,t)Vxx{x,t) + C3BVxxxx{x,t)]}, 
R T 

(3) 

where ci, C2 and C3 are constants to be determined. Integrating by parts and 
taking vt\R = Vx\r = Vx\t = ^ lead to 



J{v) = J dx J dt{v{x,t)[-ciVx{x,t)vt{x,t) - -C2Avl{x,t) + C3Bvlx{x,t)]}. 

R T 

(4) 

The unknown constants Cj (« = 1, 2, 3) can be determined by taking the 
variation of the functional (4) to make it optimal. Taking the variation of this 
functional, integrating each term by parts and make the variation optimum give 
the following relation 

2ciVxt{x,t) + 3c2Avx{x,t)vxx{x,t) + 2c3Bvxxxx{x,t) = 0. (5) 

As this equation must be equal to equation (2), the unknown constants are 
given as 

ci = 1/2, C2 = 1/3 and C3 = 1/2. (6) 

Therefore, the functional given by (6) gives the Lagrangian of the regular 
KdV equation as 

L{vt, Vx,Vxx) = ~Vx{x,t)vt{x,t) - ^Avl{x,t) + ^Bvlx{x,t). (7) 
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Similar to this form, the Lagrangian of the time-fractional version of the 
KdV equation can be written in the form 



F{oD^v, v^,v^^) - ~^[oD^v{x,t)]v^{x,t)-^Avl{x,t) + ^Bvl^{x,t), 

< a < 1, (8) 



where the fractional derivative is represented, using the left Riemann-Liouville 
fractional derivative definition as [18, 19] 



''^'^^^^^Tiihaj^t^^i dr{t-rr-'^-'f{r)],k-l<a<k,t€[a,b]. (9) 
The functional of the time-FKdV equation can be represented in the form 
J{v) = j dxj dt F{oDy, v^,v^^), (10) 

R T 

where the time-fractional Lagrangian F{oD"v, Vx,Vxx) is defined by (8). 

Following Agrawal's method [3, 4], the variation of functional (10) with 
respect to v{x,t) leads to 

SJiv) = Jdxjdt {^JoD?v + ^JVx + ^JVxx}. (11) 

B. T 

The formula for fractional integration by parts reads [3, 18, 19] 

dtf{t)aDfg{t)^ f dtg{t)tD^f{t), f{t),g{t)e[a,b]. (12) 

where t-D^, the right Riemann-Liouville fractional derivative, is defined by [18, 
19] 



^^"•^^^^ = Y(k^^^l dr{T-tf—'f{T)], k-l<a<k,te [a,b]. (13) 
Integrating the right-hand side of (11) by parts using formula (12) leads to 



R T 

where it is assumed that 5v\t = 5v\r = 5vx\r = 0. 

Optimizing this variation of the functional J{v), i. e; 5J{v) = 0, gives the 
Euler-Lagrange equation for the time-FKdV equation in the form 
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) = 0. 



(15) 



Substituting the Lagrangian of the time-FKdV equation (8) into this Euler- 
Lagrange formula (15) gives 



Substituting for the potential function, Vx{x,t) = u{x,t), gives the time- 
FKdV equation for the state function u{x, t) in the form 



where the fractional derivatives oD" and t-D^^ are, respectively the left and 
right Ricmann-Liouvillc fractional derivatives and are defined by (9) and (13). 
The time-FKdV equation represented in (17) can be rewritten by the formula 



where the fractional operator q'D^ is called Riesz fractional derivative and can 
be represented by [4, 18, 19] 



The nonlinear fractional differential equations have been solved using dif- 
ferent techniques [18-23]. In this paper, a variational-iteration method (VIM) 
[24, 25] has been used to solve the time-FKdV equation that formulated using 
Euler-Lagrange variational technique. 

3 Variational-iteration Method 

Variational-iteration method (VIM) [25-27] has been used successfully to solve 
different types of integer nonlinear differential equations [31-35]. Also, VIM 
is used to solve linear and nonlinear fractional differential equations [25, 36- 
38]. This VIM has been used in this paper to solve the formulated time-FKdV 
equation. 

A general Lagrange multiplier method is constructed to solve non-linear 
problems, which was first proposed to solve problems in quantum mechanics 







k-1 < a < k, t e [a,b]. 



(19) 
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[25]. The VIM is a modification of tliis Lagrange multiplier method [26, 27]. The 
basic features of the VIM are as follows. The solution of a linear mathematical 
problem or the initial (boundary) condition of the nonlinear problem is used as 
initial approximation or trail function. A more highly precise approximation 
can be obtained using iteration correction functional. Considering a nonlinear 

partial differential equation consists of a linear part LU{x,t), nonlinear part 

NU{x,t) and a free term f{x,t) represented as 

LU{x,t) + NU{x,t) = f{x,t), (20) 

where L is the linear operator and A'^ is the nonlinear operator. According 
to the VIM, the (n + l)th approximation solution of (20) can be given by the 
iteration correction functional as [24, 25] 



Un+l{x,t)^Unix,t)+ / dTX{T)[LUn{x,T)+NUn{x,T)- f{x,T)],n>0, 

Jo 

(21) 

where A(t) is a Lagrangian multiplier and Un{x,T) is considered as a restricted 

variation function, i. e; 5Un{x,T) = 0. Extreme the variation of the correction 
functional (21) leads to the Lagrangian multiplier A(t). The initial iteration 
can be used as the solution of the linear part of (20) or the initial value U{x, 0). 
As n tends to infinity, the iteration leads to the exact solution of (20), i. e; 



U{xJ) = lim Un{x,t). (22) 

For linear problems, the exact solution can be given using this method in 
only one step where its Lagrangian multiplier can be exactly identified. 



4 Time-fractional KdV equation Solution 

The time-FKdV equation represented by (18) can be solved using the VIM by 

the iteration correction functional (21) as follows: 

Affecting from left by the fractional operator on (18) leads to 



^u{x,t) = oD,"-\{x,t)\t=o ^^^_^^ 

- ^Dt^-"[A u{x,t)^u{x,t) + B ■^u{x,t)], 
< a < 1, i e [0,To], (23) 

where the following fractional derivative property is used [18, 19] 
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k-l < p <k. (24) 

As a < 1, the Riesz fractional derivative (f_Dj"~^ is considered as Riesz 
fractional integral that is defined by [18, 19] 



PtV{t) = lioirm + thVm = IyIo) i dr\t-rr-'f{T), a > 0. (25) 

where qI t"' f{i) and tlb"f{t) are the left and right Riemann-Liouvulle fractional 
integrals, respectively [18, 19]. 

The iterative correction functional of equation (23) is given as 



/■* d 
Un+i{x,t) = Un{x,t) + J dTA(r){— u„(x, r) 



a-2 



T(a-l) 

+ ?£>/-"[^ Un{x,T)-^Un{x,T)+B ^Un{x,T)]}, (26) 



where n > and the function Un{x,t) is considered as a restricted variation 
function, i. e; 5un{x,t) = 0. The extreme of the variation of (26) using the 
restricted variation function leads to 



SUn+l{x,t) = 6Un{x,t)+ / dTX{T) S — Un{x,T) 

Jo OT 

/■* d 

= 5Un{x,t) + X{t) SUn{x,t) - (It — X{t) 5Un{x , t) = . 

Jo OT 

This relation leads to the stationary conditions 1 + X{t) = and ^A(t) = 0, 
which leads to the Lagrangian multiplier as A(t) = — 1. 

Therefore, the correction functional (32) is given by the form 



Un+i{x,t) = 



/■* d 
in{x,t)-J dT{—U„{x,T) 

- ^I^-"Un{x,T)\r=0-, 



T(a- 1) 



d 



g3 



+ ^D^'-''[AUn{x,T)—Un{x,T) + B-^Un{x,T)]}, (27) 



dx 
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where n > 0. 

In Physics, if t denotes the time- variable, the right Riemann-Liouville frac- 
tional derivative is interpreted as a future state of the process. For this reason, 
the right-derivative is usiially neglected in applications, when the present state 
of the process does not depend on the results of the future development [3]. 
Therefore, the right-derivative is used equal to zero in the following calcula- 
tions. 

The zero order correction of the solution can be taken as the initial value of 
the state variable, which is taken in this case as 

uo{x,t) = u{x,0) = Aosech.'^{cx). (28) 
where and c are constants. 

Substituting this zero order approximation into (27) and using the definition 
of the fractional derivative (19) lead to the first order approximation as 



ui{x,t) = j4o sech^(ca;) 

+2Aoc sinh(ca;) sech^(ca;)[4c^B 

+(^o^-12c^g)sech^(cx)] ^/" (29) 

1 (a -I- Ij 

Substituting this equation into (27), using the definition (19) and the Maple 
package lead to the second order approximation in the form 



U2{x,t) = Aosech^(cx) 

+Aqc sinh(cx) sec h^(cx) 

* [Ac^B -h (^0-4 - 12c^B) sec h^ (cx)] - * 



T(a + 1) 
+Aoc'^ sec h^(cx) 

*[16c'*B2 + 8c^B{5AoA - 63c^B) scch^(cx) 
+ i3AlA'^ - UeAoc^AB + 168c^B^) sech'*(cx) 

— {AlA^ - 42Aoc^AB + SGOc'B^) sech«(cx)]^,^^^^ 

-t-AgC'' sinh(cx) sech^(ca;) 

^iAc^Blic^B + 3{AoA - Uc^B) sech^(cx)] 

+2iAlA'^ - 32Aoc^AB + 2mc^B'^) sech'^(cx) 

r(2a 1) t^a 



^(A^A2 - 24.Ao^AB + UAc'^B^) sech^(cx)} 



[r(a + l)]2r(3a + l)' 

The higher order approximations can be calculated using the Maple or the 
Mathematica package to the appropriate order where the infinite approximation 
leads to the exact solution. 



8 



5 Results and calculations 



The calculations are carried out for the solution of the time-FKdV equation 
using the VIM for different values of the equation parameters of nonlinearity 
(A) and dispersion (B). The initial value of the solution for all eases is taken as 
Aq scch(c.T) where the amplitude Aq is taken equal to unity and the constant (c) 
is used with different values as B/A, B/{4A), B/{8A) or B/{12A). The solution 
is calculated for different values of the fractional order a {a — 3/4, a = 1/2, 
a ^ 1/3 and a = 1/4). 

The 3-dimensional representation of the solution of the time-FKdV equation 
with space x and time t for constant c ~ B / (4^4) and A — \, B = \ for different 
values of the fractional order (a) is given in Fig (1). While in Fig (2), the 3- 
dimensional representation of the solution of time-FKdV equation with A — 1, 
B = 1 is given for the fractional order a = 1/2 for different values of the constant 
c = B/A, c = BI{AA), c = -B/(8A) and c = B/(12A). 

Figure (3) represents the solution of the time-FKdV equation for fractional 
order a = 1/2 with c = B/{4A) using {A = I, B = 1), {A = l, B = 3), {A = 6, 
B = 1) and {A = 6, B = 3). 

In Fig (4), the solution of the time-FKdV equation u{x,t), for {A = 1, 
B = 1) with c = B/{4:A), is calculated for a = 1/2. These calculations are 
represented as 2-dimensional figure against the space x for different values of 
time t. 
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Figure Captions 

Fig. 1: The distribution function u{x, t) with A= 1, B = 1, initial condition 
uq = sech^(cx) and c = B/(4A) for different values of the fractional order (a). 

Fig. 2: The distribution function u{x, t) with A = 1, B = 1, initial condition 
Uo = sech-^(ca;) and fractional order a = 1/2 for different values of the constant 
(c). 

Fig. 3: The distribution function u{x,t) with fractional order a = 1/2, 
initial condition uq = sech^(ca;) and c = B/{4:A) for different values of A and 
B. 

Fig. 4: The distribution function u{x,t) as a function of space (x) for 
different values of time t with A = 1, B = 1, initial condition uq = sech^(ca;), 
c = B/{4A) and a = 1/2 . 
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